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Using quasiparticle models and imposing thermodynamic consistency, lattice data for the 
equation of state of deconfined QCD can be mapped to finite chemical potential. We con- 
sider a refinement of existing simple massive quasiparticle models using the non-local hard- 
thermal-loop (HTL) propagators, and certain NLO corrections thereof, to obtain the ther- 
modynamic potential as a function of temperature and chemical potential. At small chemical 
potential we find that the results for the slope of constant pressure from our main model for 



' 2 massless quark flavors is in good agreement with recent lattice data for 2+1 flavors while it 



deviates somewhat from lattice data for 2 flavors from another group. For zero temperature, 
we obtain an estimate for the critical chemical potential which is close to that obtained from 
simpler quasiparticle models. 



I. INTRODUCTION 

■^J- , Asymptotic freedom suggests that at sufficiently high temperature and/or quark chemical po- 

tential QCD is deconfined, i.e. can be described in terms of the fundamental quark and gluonic 
degrees of freedom 0|. High temperature and small chemical potential are of relevance to the quest 
for the quark-gluon plasma in heavy-ion experiments; sufficiently high chemical potential (at com- 
paratively low temperature) may be reached in the cores of compact stars. In the latter case, novel 
color superconducting phases may occur [jj, which should however have only comparatively small 
effects on the equation of state. 

The thermodynamic pressure of QCD at high temperature T <C Aqcd has been determined in 

qj \ perturbation theory to order aj 2 in QCD @,|| and even to a 3 log a in a recent heroic effort 0. 

Strict perturbation theory, however, shows extremely poor convergence for any temperature of prac- 
tical interest so that further resummations appear to be necessary. In recent years, resummations 
which are based on the hard-thermal-loop (HTL) effective action [f| have been proposed, alterna- 
tively in form of so-called HTL perturbation theory 0, @] or based upon the 2-loop ^-derivable 
approximation 0,0,0] (see also Peshier 01). The latter approach, which assumes weakly inter- 
acting quasiparticles as determined by the HTL propagators and NLO corrections thereof, leads to 
results which agree remarkably well with lattice data down to about 3T C if standard 2-loop running 
a s is adopted, 1 and it has been recently shown that these results are in fact consistent with those 
obtained in high-order perturbation theory if the latter is organized through effective (dimension- 
ally reduced) field theory according to and effective-field-theory parameters are kept without 
further expanding in powers of the coupling [ijj], as already advocated in Ref. 0|. 

The behavior of the thermodynamic potentials closer to the transition temperature needs to be 
studied by fully non-perturbative means. For zero quark chemical potential \x = lattice QCD 
calculations are seemingly up to the task of determining thermodynamic quantities for a quark- 
gluon plasma very accura tely by no w ll.4llf|. A lthough recently there has been important progress 

13 El El SMS 



also for non- vanishing \i [TJ, LL$, 12TJ, [2lJ, l22j , the case of large /i> T (relevant for the e quat ion 
of state for cold dense matter, which is of importance in astrophysical situations 0, 24], 2r1. l2fj|) is 



1 This is in fact not the case for 2-loop HTL perturbation theory Q, but it appears that the source of the difficulty has 
to do with the necessity of thermal counter-terms and incompletely compensated hard contributions, see Ref. 
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currently beyond the reach of lattice calculations. 

As a remedy for this situation, Peshier et al. 0, 13] proposed a method which can be used 
to map the available lattice data for /i = to finite pL and small temperatures by describing the 
interacting plasma as a system of massive quasiparticles (QPs); section 2 gives a short review of 
this technique. While in their simple QP model the (thermal) masses of the quarks and gluons are 
approximated by the asymptotic limit of the hard-thermal-loop (HTL) self- energies, we consider 
models based on the HTL-resummed entropy 0, 0, O], which include more of the physically 
important plasmon effect than in the simple QP model, and (NLO) extensions thereof, which include 
the full plasmon effect. In perturbation theory it is essentially this effect that spoils convergence, 
so a priori it was not granted that a model including this effect fully would give physically sensible 
results. As it turns out, however, the results based on our models (which are presented in section 
3 and 4, respectively) turn out to be well-behaved and indeed lie not too far from the simple QP 
model. 

In section 5 we give the 2-parameter fits of the various models to the lattice data for 2 flavors 
and compare the model results for the pressure at small temperatures and large chemical potential 
in section 6. Furthermore, we compare our results for the isobar emerging at 2 T c to recent lattice 
results and give our conclusions in section 7. 



II. MAPPING LATTICE DATA TO FINITE CHEMICAL POTENTIAL USING QP 

MODELS 

Assuming an SU(iV) plasma of gluons and Nf light quarks in thermodynamic equilibrium can 
be described as a weakly interacti ng g as of massive quasiparticles with residual interaction B, the 
pressure of the system is given by [23j 

p {T, M) = w ' m ^ ~ jB ( m f (!) 

i=g,q 

where the sum runs over gluons (g), quarks and anti-quarks (q) with respective chemical potential 
0, ±/i; pi are the model dependent QP pressures and m, are the QP masses. The latter are functions 
of an effective strong coupling G 2 (T, fj,), which appears only within TOj(G 2 ,T, //) in a predefined 
form. 

Using the stationarity of the thermodynamic potential under variation of the self-energies and 
Maxwell's relations j^l, one obtains a partial differential equation for G 2 , 

dG 2 dG 2 

flT 5T +a ^ = 6 ' (2) 

where ax, and b are coefficients that are given by integrals depending on T, p and G 2 . Given 
a valid boundary condition, a solution for G 2 (T,/j,) is found by solving the above flow equation by 
the method of characteristics; once G 2 is thus known in the T, \i plane, the QP pressure is fixed 
completely. The residual interaction B is then given by the integral 

where Bq is an integration constant that has to be fixed by lattice data (usually by requiring 
P(T c ,n = 0) = Piattice{T c )) ■ Motivated by the fact that at \i = and T> T c the coupling should 



2 Here and in the following T c is always to be understood as T c | M=i 



3 



behave as predicted by the perturbative QCD beta-function, Ref. |23j used the ansatz 

487T 2 

(11JV — 2iV/)ln^iA 



° 2 ( T < °) = o at m__ T+T a \ ■ ( 4 ) 



The parameters A and T s are determined by fitting the entropy of the model (which is independent 
of B) to available lattice data at \i = 0. Using J3J) as boundary condition for Q, the above 
procedure allows one to map the lattice pressure from \x = to the whole T, //-plane. 



III. SIMPLE AND HTL QP MODEL 



A simple ansatz for the QP pressure is obtained by taking pi to be the pressure of a free gas of 
massive particles with the masses defined as the asymptotic (i.e. large momentum) gluonic 

and fermionic self-energies, respectively 



Mi 



(2N + N f ) 
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G 2 T 2 



f 4tt 2 
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(5) 
(6) 



Comparing with perturbation theory one finds that while the Stefan-Boltzmann and leading- 
order interaction terms are correctly reproduced, only l/(4\/2) of the NLO term of the interaction 
pressure (the plasmon effect ~ G 3 T 4 ) is present (lfil |. 

To include more of this physically relevant effect, one can make use of the HTL-resummed 
entropy 0, 0, O], which takes into account the momentum dependence of the QP excitations as 
well as Landau-damping effects. The HTL model ansatz for the QP pressure for gluons (p g ) and 
quarks p q then reads 



Pg 



Pq 



d 3 k 



n(io) 



21m In -uj 2 + k 2 + n T - 2Imn T ReD T 



Im In k 2 + n L + Imn L ReD L 
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d 3 k r°° (to 
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(7) 



(8) 



where d g = 2(N 2 — l),d q = 2NNf for gluons and quarks/anti-quarks, respectively; n(uj) and f±(uj) 
are the bosonic and fermionic distribution functions. Dt,l, A-t are the HTL propagators with LTt,l 
and S± the corresponding self- energies, 
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(9) 
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(11) 



which we have written in terms of the HTL Debye mass squared m 2 D = 27/1^ and the HTL zero- 
momentum quark mass squared M 2 = M^/2. 
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The HTL QP model includes 1/4 of the full plasmon effect at [i = 0, which is a factor of v2 
more than in the simple QP model which only uses the asymptotic HTL masses. The remainder of 
the plasmon effect is in fact entirely due to NLO corrections to these asymptotic masses at order 
G 3 T 2 UM\H 



IV. NLO MODELS 

The complete (momentum-dependent) NLO corrections to the asymptotic masses of quarks and 
gluons have not been calculated yet, but as far as the plasmon effect is concerned, these corrections 
contribute in the averaged form p, llfj] 

- 9 f dkkn'(k)Re5U T (Lo = k) 

Sm » = J TWT = NTmD (12) 

and similarly 

6M ~ = Jdkk(fUk) + fL{k)) = 2, • (13) 

For the values of the coupling G considered here, these corrections are so large that they give 
tachyonic masses when treated strictly perturbatively. In Ref. 0] it has been proposed to incor- 
porate these corrections through a quadratic gap equation which works well as an approximation 
in the exactly solvable scalar 0(N — > oo)-model, where strict perturbation theory would lead to 
identical difficulties. However, for the fermionic asymptotic masses, in order to have the correct 
scaling of Casimir factors in the exactly solvable large- iVj limit of QCD |29|, a corresponding gap 
equation has to remain linear in the fermionic mass squared. Choosing the correction therein to be 
determined by the solution to the gluonic gap leads to 0, 

_ 2 . 2 G 2 NT _ 

m oo = m oo ~ m oo (14) 

Ml = Ml- G ^ N2 ~ ^K moo , (15) 

where ml and Ml are the leading-order gluonic and fermionic asymptotic masses as given in Q 
and (0. This in fact avoids tachyonic masses for the fermions as long as Nf < 3. 3 

Finally, since the averaged quantities ml and Ml are the effective masses at hard momenta 
only, a cutoff scale A = \j2TTTrriE>c\ is introduced that separates soft from hard momenta. The QP 
pressure for this model, which in the following will be referred to as NLA-model, then separates into 
a soft and a hard component for both gluons and fermions. The soft parts are given by expressions 
like and iJSJ), but with A as upper limit for the momentum integration. For the hard parts, 
the momentum integrations run from A to oo and the mass pre-factors in the HTL self-energies 
I19I10I11II are replaced by their asymptotic counterparts, m 2 D — ► 2ml an d M 2 — > \Ml. 

The single free parameter c\ in A can be varied around 1 to obtain an idea of the "theoretical 
error" of the model. In the following we will consider the range ca = \ to c\ = 4; note that ca = oo 
corresponds to the HTL-model (which has been used as cross-check) since all hard corrections are 



3 For Nf = 3 the solutions to the above approximate gap equations happen to coincide with those obtained in the 
original version of two independent quadratic gap equations of Ref. 10]. For the case Nf — 2 considered in this 
paper, the differences are fairly small. For Nf > 3, however, the necessity to avoid tachyonic masses would restrict 
the range of permissible coupling strength G. 
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FIG. 1: Entropy data generated from Ref. vs. 
fitted model entropy. The band was generated from 
NLA model for ca between 4 and 1; the HTL and 
simple QP models lie at the lower boundary. 
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FIG. 2: Effective coupling: NLA model results for 
ca = 4,1 and 1/4 (full lines from bottom to top), 2- 
loop perturbative coupling in MS (gray band) and 
result from Ref. [3] (dotted line). 



ignored. On the other hand, ca = would assume that lfT4"|) and l|T5jl represent good approximations 
for the NLO corrections to the spectral properties of soft excitations. However, the few exis ting 
results, in particular on NLO corrections to the Debye mass and the plasma frequency [22], 
appear to be rather different so that it seems safer to leave the soft sector unchanged by keeping a 
finite ca- 



V. ENTROPY AND PRESSURE FOR \i = 

To obtain the input parameters T* and A we fitted the entropy expressions from the models 
under consideration to lattice data |2j]] for Nf = 2, with an estimated continuum extrapolation as 
used in Ref. Q. We find 





HTL 


c A = 4 


C A = 1 


ca = 1/4 


T s /T c 


-0.89 


-0.89 


-0.84 


-0.61 


A 


19.4 


18.64 


11.43 


3.43 



where it can be seen that the NLA model with ca = 4 is very close to the HTL model. The fits to 
the entropy data (shown in Fig.^) all lie in a narrow band for all the models considered. The fitted 
effective coupling G 2 for the various models is shown in figure [2j for comparison, also the 2-loop 
perturbative running coupling in MS is shown, where the renormalization scale is varied between 
7rT and 4ttT and following Ref. we have chosen T c = 0.49Aj^g. As can be seen from the plot, 
the results for the effective coupling are well within the range of the 2-loop perturbative running 
coupling (for the case ca = 1/4 and renormalization scale irT the results even seem to be identical, 
which is, however, probably only a coincidence). We also show the result for the coupling obtained 
in the semiclassical approach of Ref. j3| 4 . 

In general one can see that the effective coupling becomes bigger when ca gets smaller; this is 
because the hard masses (for equal values of the coupling) are smaller than the soft masses which 



4 For finite fi and constant temperatures, however, the coupling obtained in [^U rises while our results indicate a 
decrease of the coupling (which is consistent with the standard QCD running coupling with renormalization scale 
proportional to \JT 2 + (/j,/tv) 2 ). 
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B(T,/z = 0)/T 4 




FIG. 3: Residual interaction B: NLA model results 
between ca = 4, ca = 1 (full line) and ca = 1/4 
(lower boundary). 
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FIG. 4: Pressure data from Ref.[27] vs. model pres- 
sure. The band shows the NLA model results for ca 
between 4 and 1/4. 



makes the entropy increase when the hard parts become more important. Accordingly, the coupling 
has to rise in order for the entropy to match the data (therefore, for the extreme case ca = one 
finds huge values of the effective coupling constant). 

Once the effective coupling is known for \i = one can proceed to evaluate the QP pressure 
Pi and the residual interaction B from Eq. (JHJ, shown in Fig. in this figure, the HTL results 
(not shown) would lie marginally above the upper boundary of the band (cf. |35l | for a compar- 
ison between simple QP and HTL model). A comparison between the full pressure £Q| and the 
lattice data is shown in figureQJ where the integration constant Bq was set so that P(T C ) = 0.536(1): 





HTL 


c A = 4 


C A = 1 


ca = 1/4 




0.82 


0.73 


0.6 


0.47 



VI. FINITE CHEMICAL POTENTIAL 

A calculation and subsequent evaluation of the coefficients a^, ax and b in shows that b, in 
contrast to a M and ay, is strongly model-dependent. Solving the flow equation (J2J and integrating 
B along the characteristics using Eq. © one obtains pressure and effective coupling in the whole 
T, \i plane. Here we focus on the pressure but the calculation of the other thermodynamic quantities 
is straightforward once the effective coupling and the pressure have been obtained. 

Concerning the shape of the characteristics, it has already been noticed in jijj that the crossing 
of characteristics in the simple QP model [23J does not occur in the HTL-model; this feature is 
preserved in the NLA models. 



A. Lines of constant pressure 

It is straightforward to extract from our data the line where P(T,fi) = P(T c ,0); for Nf = 2, 
the results for the simple and HTL QP model as well as for the NLA models for ca > 1 are rather 
close to each other. For NLA models with ca < 1, however, the slope of the constant pressure line 
(until n/T c = 1) deviates rapidly from the values found for ca > 1: 
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FIG. 5: Lines of constant pressure: NLA models 
for Nf — 2 from c\ — 4, ca = 1 and ca = 1/4 
(dotted lines from lowest to highest); lattice data 
for N f = 2 + 1 [ioj (light gray band) and Nf = 2 
[Ls| (long dashed- lines). 
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FIG. 6: Difference £P scaled with T 4 for m/^ c = 1. 
Shown are HTL model (triangles) and NLA models 
(boxes). See text for details. 





HTL 


c A = 4 


c A = 1 


c A = 1/4 


rp dT 


-0.06818(8) 


-0.06810(6) 


-0.06329(34) 


-0.041(9) 



We have also compared these results with recent lattice data for 2+1 flavors [20] in figure El The 
variations of the lattice data in this figure result from fitting the data with a 2nd order (upper limit) 
and 4th order polynomial (lower limit) as used in j^|. The comparison shows that our 2-flavor 
results obtained from simple, HTL and NLA QP models with ca > 1 (indicated in figure 03 by the 
two lower lying pointed lines) are in very good agreement with the lattice data whereas for NLA 
models with ca < 1 the slope is much natter. Also shown are the lattice results for the constant 
pressure line from Q for 2 flavors, with a slope of T c |^ = -0.107(22), which is significantly 
steeper. 



B. Susceptibilities 

It been noticed in |2q| that the quantity 

AP P(T,fj,)-P(T,0) 



AP SB Psb(T,ii)-Psb(T,0) 



(16) 



is essentially /i-independent; we recover a similar scaling behavior for our models. As expected, for 
small chemical potential this curve is very close to the quark-number susceptibilities x(T)/xo(T) 
obtained at fi = 0. In fact, we found that also for larger \i the pressure can be well approximated 
by 

P(T,/i) - P(T,0) = ^(/x 2 + ^2). (17) 

To quantify this assertion we consider the subtracted quantity SP = P(T, /i) — P(T, 0) — ^y^(/i 2 + 

4 

2^± T -± )) which should be zero, if the pressure at finite chemical potential was fully determined by 
the susceptibilities at fi = according to (fT7|l : the results (scaled with T 4 ) are shown in figureElfor 
n/T c ~ 1. As can be seen, dropping the fi 4 term in Eq. lfT7)l (open symbols) deteriorates the result 



X/Xo 



P/Ps 



0.8 



. 4 



0.2 - 




n — 



1.5 



3 

T/T c 



FIG. 7: Susceptibilities from NLA models (light 
gray band), from Nf — 2 lattice data (trian- 
gles) and "scaling curve" for 2+1 flavors from [2(il | 
(boxes). 
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FIG. 8: The pressure at vanishing T: Shown are per- 
turbative results (light gray band) and NLA models 
for ca from 1/4 to 4 (dark band). 
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considerably. Trying to determine the higher-order susceptibility, % — -g-. 
as 

4! 



which we parameterized 



X 



X(T) 



2T 2 x{T f 



we find that - within our numerical errors - the result for x(T) is consistent with y/2ir = 4.44 . . . 
for T > 1.5T C (which was our original assumption in Eq. ) ; therefore, our results are close to 
but about 10 percent lower than the values found in Ref. [22]. 

Clearly, 5P gets smaller with increasing temperature, but also at relatively low temperatures 
and high chemical potential, 5P is remarkably small (e.g., for the HTL model for T = 1.5T C and 
fi/T c ~ 4 we find a 5P/T 4 smaller than 10 percent). 

In Fig. we compare the results for the susceptibilities of our NLA models with the lattice data 
for 2 flavors from Ref. and the "scaling curve" in 0|, which provides an approximation to the 
2+1 flavor quark-number susceptibility at \x = 0. We find that all these results agree very well, 
with the agreement getting better at higher temperatures. 

It should be noted however that none of the lattice data we are using has yet been rigorously 
extrapolated to the continuum limit, so that both the data and the extrapolation by means of QP 
models are still likely to change somewhat when this will be done eventually. 



C. Very small temperatures 

Extending the lines of constant pressure from our models to very small temperatures we obtain 
a crude estimate of the phase transition line in this region of phase space. Denoting with [a c the 
chemical potential where (for vanishing temperature) the pressure equals the lattice pressure at 
Li = 0, P(0,fj, c ) = P(T C , 0), we find for our models (assuming T c = 172 MeV) 





HTL 


c A = 4 


c A = 1 


c A = 1/4 




533 MeV 


536 MeV 


558 MeV 


584 MeV 


Ho 


509 MeV 


511 MeV 


537 MeV 


567 MeV 



Here we have also given the values (Xq where the pressure vanishes, which may be taken as a definite 
lower bound for the critical chemical potential within the respective models. 
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In general, these results are in agreement with the estimates for \i c from |2J, |26j; for the QP 
models considered, the lowest and highest results for \x c are obtained for the HTL model and the 
NLA model with c\ = 1/4, respectively, while the /x c of the simple QP model lies between the NLA 
ca = 4 and c\ = 1 model values. However, even our lowest result for \i c turns out to exceed the 
value for the critical chemical potential expected in Ref. E^. 

A comparison of the perturbative pressure at vanishing temperature with the results for our 
models at T/T c ~ 0.01 is shown in figure |H1 We plotted the perturbative results converted to MS 
with the standard 2-loop running coupling, Aj^g = T c /0.49 [33] and renormalization scale varied 
from n to 3/j,; a comparison of this coupling and our results is shown in figure El 



D. Equation of state for cold dense matter 



By calculating the number density at T/T c ~ 0.01 and using our results for the pressure we 
obtain an equation of state for cold deconfined matter. As is the case for the simple QP model Ei" 
we find that the energy density S is well fitted by the linear relation 



8{p) =AB + ap 



with 





HTL 


c A = 4 


c A = 1 


c A = 1/4 


4B/r c 4 


11.1(8) 


12.3(8) 


14.7(9) 


19.2(1.6) 


a 


3.23(5) 


3.22(4) 


3.22(4) 


3.17(4) 



A value of a ~ 3.2 for Nf = 2 has also been found in the simple quasiparticle model in [24] and 
thus seems to be model independent, in contrast to the bag constant .B 1 / 4 , which varies between 
314 and 360 MeV. By using the Tolman-Oppenheimer-Volkov equations and the equation of state 
one can determine the mass-radius relations of non-rotating quark-stars j^]; choosing a = 3.2 and 
taking the bag constant values from above, we find stable star configurations with radii ranging 
from 3.6 to 4.9 km and masses from 0.7-0.8 solar masses for NLA, ca = 1/4 and HTL models, 
respectively. However, it should be kept in mind that the outermost layers of such a quark star 
are metastable with respect to hadronic matter and therefore the details of the star structure will 
depend sensitively on the hadronic equation of state I2T 



VII. CONCLUSIONS AND OUTLOOK 

We considered an improvement of simple quasiparticle models (^J by using the full HTL ap- 
proximation and certain NLO corrections thereof along the lines of Ref. |9j, LLjj, lll| and investigated 
the respective predictions at finite chemical potential when these models are matched to Nf = 2 
lattice data at fx = 0. 

We found that the slope of constant pressure for those models with ca > 1 agrees very well 
with recent lattice data for 2+1 flavors [i^, EH], while the lattice data for Nf = 2 indicate 
significantly steeper slopes. 

We also found that for our models the pressure for non- vanishing \x scales with the quark number 
susceptibilities at fi = 0, provided that the // 4 terms (corresponding to a higher order susceptibility) 
are not dropped but added with the same weight as in the ideal gas limit. We therefore interpreted 
the "scaling curve" from 2+1 flavor lattice data [^J as susceptibilities at \i = and found that a 
comparison with Nf = 2 susceptibilities from another lattice group as well as those obtained 
from our models all agree very well. 
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FIG. 9: Effective Coupling from NLA models with 
ca = 4 to 1/4 (dark band) and perturbative 2-loop 
running coupling (light gray band). 



We finally extended our results to the small temperature, high density regime and obtained an 
estimate for the critical de nsity at T = which is consistent with earlier results, though in excess 
of the expectations of Ref. jiM Issl]. Furthermore, we obtained an equation of state for cold, dense 
matter, which allows for pure quarks stars with masses of ~ O.8M and radii of less than 5 km, 
similar to the results obtained by 
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